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ABSTRACT 

E. V. Pitjeva, by processing more than 400,000 planetary observations of var- 
ious types with the dynamical models of the EPM2006 ephemerides, recently 
estimated a correction to the canonical Newtonian-Einsteinian Venus' perihelion 
precession of —0.0004 ±0.0001 arcseconds per century. The prediction of general 
relativity for the Lense-Thirring precession of the perihelion of Venus is —0.0003 
arcseconds per century. It turns out that neither other mismodelld/unmodelled 
standard Newtonian/Einsteinian effects nor exotic ones, postulated to, e.g., ex- 
plain the Pioneer anomaly, may have caused the determined extra-precession of 
the Venus orbit which, thus, can be reasonably attributed to the gravitomag- 
netic field of the Sun, not modelled in the routines of the EPM2006 ephemerides. 
However, it must be noted that the quoted error is the formal, statistical one; the 
realistic uncertainty might be larger. Future improvements of the inner planets' 
ephemerides, with the inclusion of the Messenger and Venus-Express tracking 
data, should further improve the accuracy and the consistency of such a test of 
general relativity which would also benefit of the independent estimation of the 
extra-precessions of the perihelia (and the nodes) by other teams of astronomers. 

Subject headings: Experimental tests of gravitational theories; Celestial mechanics; 
Orbit determination and improvement; Ephemerides, almanacs, and calendars 
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In the weak-field and slow motion approximation the Einstein field equations of general 



relativity get linearized resembling the Maxw ellian equations o : 



m 



2001 



electromagnetism. As 



20071 ); it is induced by 



a consequence, a gravitomagnetic field arises flMashhoor 
the off-diagonal components goi,i = 1,2,3 of the space-time metric tensor related to the 
mass-energy currents of the source of the gravitational field. It affects orbiting test particles, 
precessing gyroscopes, moving clocks and at oms and propagating electromagnetic waves 



(JRuggiero and Tartaglia 



2002 



S chafer 



20041 ). The most famous gravitomagneti c effects 



1959 



Schifi 



19601 ). whose 



are, perhaps, the precession of the axis of a gyroscope flPugh 
detection in t he gravitational fie ld of the rotating Earth is the goal of t he space-based GP-B 
exper iment!^ (jEveritt et all Il974l ) , and the Lense-Thirringp precessions (ILense and Thirring 
19181 ) of the orbit of a test particle for wh ich some disputed satellit e -based tests in th e 



gravitational fields of the sp inning EarthJ 



2006a 



2007a 



Lucchesi 



20051 ) and Mars fllorio 



Ciufolini and Pavlis 



2006b ; 



Krogh 



2004 



2007 



Ioriol 



2005; 



Iorio 



2005a 



2008af ) have been 



reported. 

We focus on the detection of the solar gravitomagnetic field through the Lense-Thirring 
planetary precessions of the longitudes of periheliarj w = to + cos i f2 

dm AGScosi 

where G is the Newtonian gravitational constant, S is the proper angular momentum of the 
Sun, c is the speed of light in vacuum, a and e are the semimajor axis and the eccentricity, 



1 See on the WEB http://einstein.stanford.edu/ 



20071 ). it 



2 According to an interesting historical analysis recently performed in (jPfister 
would be more correct to speak about an Einstein-Thirring-Lense effect. 

3 Here uj is the argument of pericentre, reckoned from the line of the nodes, i is the 
inclination of the orbital plane to the equator of the central rotating mass and Q is the 



longitude of the ascending node. 
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respe ctively, of the planet's orbit. It may be interesting to know that in (IHaas and Ross 



19751 ) it was proposed to measure the solar gravitomagnetic field through the Schiff effect 



with a drag- free gyroscope orbiting the Sun in a polar orbit. 

The impact of the Sun's rotation on th e Mercury's long itude of perihelion was calculated 



for the first time with general relativity by 



de Sitter! (119161 ) who, by assuming a homogenous 



and uniformly rotating Sun, found a secular rate of —0. 01 arc s econd s per century (" cy 1 in 



Soffell fll989h . 



Cugusi and Proverbio 



the fo llowing;) . This value is also quoted at pag. Ill of 

( 119781 ) yield —0.02 " cy^ 1 for the argument of perihelion of Mercury. Instead, recent 

determinations of the Su n's proper angular m omentunQ S & = (190.0 ± 1.5) x 10 39 kg m 2 



s 1 from helioseismology (jPijpersi Il998l . 120031 ). accurate to 0.8%, yield a precessional effect 



one order of magnitude smaller. The predicted gravitomagnetic precessions of the four 
inner planets, according to the recent value of the Sun's angular momentum, are reported 
in Table [U they are of the order of 10~ 3 — 10 -5 " cy^ 1 . Due to their extreme smallness it 
has been believed for a long time, until recently, t hat the plan etary Lense-Thirring effect 
would have been undetectable; see, e.g., p. 23 of 
showing that recent advance s in the eph e merides 



Soffell (J1989I ). A preliminary analysis 



favorable was carried out in ( Iorio 



2005a). 



i eld are making the situation more 



Pitjeva 



(j2005al ) processed more than 317,000 



planetary observations of various kinds collected from 1917 to 2003 with the dynamical 



'It could me measured also with a different approach (INi 



2008|). 



Table 1: Lense-Thirring precessions, in " cy -1 , of the longitudes of the perihelion w of the 
inner planets of the Solar System induced by the gravitomagnetic field of the Sun. The value 
Sq = (190.0 ±1.5) x 10 39 kg m 2 s _1 has been assumed for its angular momentum. 



Mercury Venus Earth 



Mars 



-0.0020 



-0.0003 



-0.0001 



-0.00003 
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force models of the EPM2004 ephemerides f|Pitjevall2005bl ) producing a global solution in 
which she estimated, among many other parameters, also corrections Am to the canonical 
Newton-Einstein perihelion precessions for all the inner planets; since the gravitomagnetic 
force was not modelled at all, contrary to the static part of the general relativistic force 
of order 0(c~ 2 ), such corrections to the usual perihelia evolutions account, in principle, 
for the Lense-Thirring effect as well, in addition to the mismodelled parts of the standard 
Newtonian/Einsteinian precessions. Thus, the estim ated correctio ns for the perihelion rates 



of Mercury, the Earth and M ars have been u sed in fjlorio 



The errors 5 (Am) released in (IPitjeva 



2007bl ) to perform a first test. 



2005al ) were slightly larger than the gravitomagnetic 



precessions whose predicted values, however, were found compatible with the estimated 
corrections. Venus was not used because of the poor data set used in the estimation of 
its extra-precession whose value, indeed, turned out too large to be due to a physically 
plausible effect amounting to +0.53 ± 0.30" cy" 1 : the Lense-Thirring prediction for the 
Venus perihelion precession was incompatible with such a result at about 2 — a level. 



Now the situation for the second planet of the Solar System has remarkab 



Pitieval (12007 



y improved 



2008|) 



allowing for a more stringent test of the Lense-Thirring effect. Indeed 
in the effort of continuously improving the planetary ephemerides, recently processed more 
than 400,000 data points (1913-2006) with the EPM2006 ephemerides which encompasses 
better dynamical models with the exception, again, of the gravitomagnetic force itself. Also 
in this case she estimated, among more than 230 par ameters, the corrections to the usual 



perihelion precessions for some pl anets 



the radiometric data of Magellan (IPitjeva 



Pitieva 



20071 ) . In the case of Venus t 



re inclusion of 



20081 ) as well allowed her to obtain_ 



Arovenus = -0.0004 ± 0.0001 " cy" 1 , 



(2) 



in which the quoted uncertainty is the formal, statistical one. By looking at Table [T] it turns 



5 Personal communication by Pitjeva to the author, June 2008. 
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out that such an extra-precession can be well accommodated by the general relativistic 
prediction for the Lense-Thirring rate of the Venus'perihelion whose existence would, thus, 
be confirmed at 25%. Somebody may object that the gravitomagnetic force should have 
been explicitly modelled and an ad-hoc parameter accounting for it should have been 
inserted in the set of parameters to be estimated. Certainly, it may be an alternative 
approach which could be implemented in future; in addition, we note that the procedure 
followed by Pitjeva may be regarded, in a certain sense, as safer for our purposes because it 
is truly model-independent and, since her goal in estimating Am was not the measurement 
of the Lense-Thirring effect, there is a priori no risk that, knowing in advance the desired 
answer, something was driven just towards the expected outcome. 

The main question to be asked is, at this point, the following one: Can the result 
of eq. ([2D be explained by other unmodelled/mismodelled canonical or non-conventional 
dynamical effects? Let us, first, examine some standard candidates like, e.g., the residual 



prec ession due to the sti 
J 2 (IPireaux and Rozelot 



imp erfect knowledge of the Sun's quadrupol e mass moment 



Pitjeva 



((2005 



20031 ) whose action was, in fact, modelled by 
by keeping it fixed to J® = 2 x 10~ 7 in the global solution in which she estimated the 
corrections to the perihelion precessions. The answer is negative since the Newtonian 



secular precession due to the Sun's o 



Dlateness 



positive. Indeed, it is (ICapderou 



2005 



Iorio 



jwhatever magnitude J 2 may have, is 



2007bJ) 



nJ 2 



W .h = o" 



3 2^2 



— 1 sin i 

a V 2 



(3) 



2(l-e 

where n = ^GM/a 3 is the Keplerian mean motion and R is the Sun's mean equatorial 
radius; the angle i between the Venus'orbit and the Sun's equator amounts tqj 3.4 deg only 
For J 2 = 2 x 10 -7 , the nominal value of the Venus' perihelion precession induced by it 



6 For an oblate body J 2 > 0. 

7 Indeed, the orbit of Venus is tilted by 3.7 deg to the mean ecliptic of J2000 
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amounts to +0.0026 " cy 1 ; by assuming an uncertainty of about 5J2 ~ 10% (IFienga et al 
20081 ). if Atbvenus was due to such a mismodelled effect it should amount to +0.0002 " 



cy 1 , which is, instead, ruled out at 6 — a level. Concerning the precession due to the solar 

_ii 



octupole mass moment J®, it is (jCapdero 



20051 ) 



wj 4 



15 (R 

16 V a 



o2\3 



;i + |e 2 ) 

(A - e 2\A 



— sin i — 2 sin 2 i-\ — 
4 5 



(4) 



For Venus it amounts to — 1.2 Jf " cy 1 ; since Jf 



Mecheri et al 



-4 x 10 9 (jRoxburgh 



2001 



20041 ). we conclude that the second even zonal harmonic of the multipolar 
expansion of the solar gravitational potential cannot be responsible for eq. ([2]) and, more 
generally, it does not represent a potentially relevant source of systematic error for the 
measurement of the Lense-Thirring planetary precessions. Similar arguments hold also 
for other potential sources of systematic errors like, e.g., the asteroid ring and the Kuiper 
Belt objects, both modelled in EPM2006: the precessions induced by them are positive. 
Indeed, a Sun-centered ring of mass m r i ng and inner and outer radius -R m m/max ~^> o induces 



a perihelion precession (jlorio 



2007ch 



W 



ring 



3 /Ga 3 (l-e 2 ) 



m 



ring 



>0. 



(5) 



IVl J^rnin^'riiax.y^'min T" -*tmax7 

According to eq. (JSJ), the precession induced by the asteroids' ring on the Venus 'perihelion 
amou nts to +0.0007 + 0.0001 " cy^ 1 by using m ring = (5 ± 1) x 1O~ 1O M flKrasinskv et al 
20021 ); the lowest value +0.0006 " cy" 1 is incompatible with eq. ((2j) at 10 — a level. In the 
case of the Kuiper B elt Objects, eq. (jSJ) yields a precession of the order of +0.00006 " cy" 1 



with m = 0.052m & (jlorio 



2007d ). Thus, we can rule out such modelled classical features 
of the Sun and the Solar System as explanations of A-cbvenus- General relativistic terms of 
order 0(c~ A ) were not modelled by Pitjeva; however, the first correction of order (9(c~ 4 ) to 



(http://ssd.jpl.nasa.gov/txt/aprx_pos_planets.pdf), while the Carrington's angle between the 



Sun's equator and the ecliptic is 7.15 deg (IBeck and Giles 



2005J). 



the perihelion precession (JDamour and Schaferlll988l ) can be safely neglected because for 
Venus it is 



W c i oc 



n{GMf 



c 4 a 2 (l 



3 2^|2 



1(T 7 "cy^ 1 . 



(6) 



Concerning possible exotic explanations, i.e. due to some modifications of the currently 



known Newton-Einstein laws o 



recently proposed extra-forces flStandish 



gravity, 



accommodate the Pioneer anomaly (INieto 



it ma y have some interest to check some of the 



20081) w hich would be able to phenomenologically 



20061 ). All of such hypothetical new forces have 
not been modelled by Pitjeva, so that if they existed in Nature they would affect Awycnus- 
A central acceleration quadratic in the rad ial component v r of the velocity of a test particlqf 



( iJaekel and Reynaud 



2005 



Standish 



20081 ) 



A = -v 2 r H, H = 6.07 x 10~ 18 m" 1 



(7) 



would induce a retrograde perihelion precession according to (jlorioll2008bl ) 

7i\ 



w 



Tinayl 



2vT 



However, eq. ([8]) predicts a precession of —0.0016 " cy" 1 for Venus, which is ruled out by 
eq. ([2]) at 12 — a level. Another possible candidate considered in (jStandishl 120081 1 is an 
acceleration linear in the radial velocity 



14 -1 



A = -\v r \lC, /C = 7.3x 10" i4 s 



(9) 



which yields a retrograde perihelion precession (jlorio 

"2e-f: 



2008bJ) 



-- 



Ky/l 



IX 



e 2 n^m) 



< 0. 



(10) 



The prediction of eq. fflOl) for Venus is —0.1 " cy 1 , clearly incompatible with eq. (T5]). 
Should one consider a central uniform acceleration with the magnitude of the Pioneer 



^The quoted numerical value of 7i allows to reproduce the Pioneer anomaly. 
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anomalous one, i.e. A = — 8.74 x 10 m s , the exotic precession induced by it 



( llorio and Giudice 



2006 



Sanders 



^Ven = A 



20061 ) on the perihelion of of Venus would be 
: -16 " cy" 1 . 



ail - e 2 



GM 



;n) 



Another non- conventiona l effect which may be considered is the precession predicted by 



Lue and Starkmanl (120031 ) in the framew ork of the PGP multidimensional braneworld 



model by Dvali, Gabadadze and Porrati (IDvali et al 
acceleration without invoking dark energy. It is 

3c 



20001 ) proposed to explain the cosmic 



^LS = =F" 



>'"0 



0(e 2 ) « +0.0005 " cy^ 1 , 



(12) 



where the plus sign is related to the self-accelerated branch, while the minus sign is for 
the standard, Friedmann-Lemaitre-Robertson- Walker (FLRW) branch; ro ~ 5 Gpc is a 
threshold characteristic of the DGP model after which gravity would experience neat 
deviations from the Newtonian-Einsteinian behavior. As can be noted, the self-accelerated 
branch is ruled out at 9 — a level by eq. ([2D, while the FLRW case is still compatible with 
eq. (T5J) (1 — a discrepancy). By the way, apart from the f act that there are theoretical 



concerns with the DGP model (see, e.g., (jlzumi et all 120071 ) and references therein), the 



existence of both the Lue-Starkman FLRW precession and the Lense-Thirring one, implying 
a total unmodelled effect of —0.0008 " cy -1 , would be ruled out by eq. ([2j) at 4 — a level. 
As a consequence, we can conclude not only that the examined exotic modifications of the 
standard laws of gravity, not modelled by Pitjeva, are not responsible for the estimated 
At^venus, but also that their existence in the inner regions of the Solar System is falsified 
by the observations. Moreover, given the magnitudes of the hypothetical effects with the 
negative sign, it is not possible that reciprocal cancelations with the positive classical 
mismodelled precessions can explain eq. (J2J). Indeed, the sum of the latter ones is +0.0004 
" cy" 1 ; the sum of, e.g., eq. ® and eq. (JTSJ (FLRW) is -0.0021 " cy" 1 , while the sum of 
eq. (jHJ) and eq. (IT2"|) (self-accelerated branch) is —0.0011 " cy -1 . 
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Thus, we conclude that the most likely explanation for eq. ([2]) is just the general 
relativistic Lense-Thirring effect. However, caution is in order in assessing the realistic 
uncertainty in such a test because, as already stated, the released error of 0.0001 " cy _1 is 
the formal, statistical one; the realistic uncertainty might be larger. By the way, we can at 
least firmly conclude that now also in the case of Venus the general relativistic predictions 
for the Lense-Thirring effect on w are compatible with the observ ational dete rminations 



2007bh . Moreover, 



for the unmodelled perihelion precessions, contrary to the case of (lloriol 
future modelling of planetary motions should take into account the relativistic effects of 
the rotation of the Sun as well. The steady improvement in the planetary ephemerides, 
which should hopefully benefit of the radiometric data from Messenger and Venus- Express 
as well, should allow for more accurate and stringent test in the near-mid future. Of 
great significance would be if also other teams of astronomers would estimate their own 
corrections to the canonical perihelion (and also node) precessions in order to enhance the 
statistical significance and robustness of this important direct test of general relativity. 
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